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ABSTRACT 

We have used /fST/NICMOS broad-band (at 1.6 /xm) and narrow-band Paa (Arost = 1.87 /im) images 
to identify star clusters and H ii regions respectively in a sample of 8 luminous infrared galaxies (LIRGs) . 
These observations have revealed the presence of a large population of super star clusters and bright 
Hii regions. A significant fraction of the Hii regions shows Ha luminosities above that of 30 Doradus, 
the prototypical giant Hii region. The excess of extremely luminous Hii regions in LIRGs has been 
confirmed by comparison with normal galaxies observed at similar spatial resolutions. Despite the large 
numbers of identified star clusters and H ii regions in LIRGs, we only find a small fraction of coincidences, 
between 4% and 30% of the total number of detected sources. Using evolutionary synthesis models we 
have reproduced the relative fractions of young Hii regions, intermediate and old star clusters observed 
in Arp 299 and the central region NGC 3256 using a Salpeter IMF and instantaneous star formation. H ii 
regions with no detected near-infrared cluster counterpart (25 — 39% of the detected sources) represent 
the youngest sites of star formation, with ages of up to approximately 5 Myr and mostly intermediate 
mass (~ 10^ M0) ionizing clusters. For these two galaxies, and within the present detection threshold we 
can only detect coincidences (4 — 10% of the detected sources) between an H 11 region and a near-infrared 
star cluster for the most massive star clusters (~ 10^ M©) during the first 7 Myr of their evolution. If 
there is significant extinction during the first million years, we may not detect the youngest star forming 
regions, and hence the observed fractions of Hii regions and coincidences will be lower limits. The 
identified near-infrared SSCs with no detectable Paa emission represent the "old" population (53 — 66% 
of the detected sources), with ages of between 7 and 20 — 40 Myr. Older clusters possibly created in 
this or previous episodes of star formation are likely to exist in these systems but cannot be identified 
with the present detection threshold. Our study demonstrates that Paa narrow-band imaging of LIRGs 
and interacting galaxies identifies the youngest sites of star formation that could be otherwise missed by 
near-infrared broad-band continuum surveys. 

Subject headings: galaxies: ISM — ISM: HII regions — galaxies: spiral — infrared radiation — 
infrared: galaxies — galaxies: interactions — galaxies: star clusters 



1. INTRODUCTION 

Infrared galaxies (Rieke & Low 1972), and in partic- 
ular nearby luminous and ultraluminous infrared galax- 
ies (LIRGs and ULIRGs, with Ljr = 10" - lO^^ Lq and 
I/iR, > 10^^ L0, respectively) have long been recognized as 
one of the best laboratories to study the process of violent 
star formation in the Local Universe. The dust-rich envi- 
ronments of LIRGs and ULIRGs are thought to be similar 
to the conditions in which star formation occurred at high 
redshift (Ivison et al. 2000 and references therein). 

It is clear now from both the theoretical and the obser- 
vational viewpoints that interactions play a major role in 
enhancing star formation throughout galaxies. In fact, the 
fraction of interactions/mergers increases from 20 — 30% 
for LiR < 101° Lq up to 70 - 95% for Ljr > lO^^ L©. The 
process of merging with accompanying super-starbursts 
appears to be an important stage in galactic evolution, 
possibly even converting spiral galaxies into ellipticals (see 
Sanders & Mirabel 1996 for a review). Moreover, it has 
been recently suggested that most interacting/merging 
systems have gone or will go through an ultraluminous 



phase, depending on how efficiently molecular gas can be 
tunneled into the nuclear regions. This may imply that 
most of ULIRG class systems may spend a large portion 
of their lifetime in the less powerful class of LIRGs (Mur- 
phy et al. 2001). 

In last decade or so, mainly with HST observations, we 
have started to get the first glimpses into the detailed na- 
ture (that is, on scales of a few tens of parsecs) of the ex- 
treme star formation processes in LIRGs and interacting 
galaxies. The discovery, both in the optical and the near- 
infrared, of a large population of massive (young) star clus- 
ters, the so-called super star clusters (SSCs), has shown 
that a significant fraction of the massive star formation is 
occurring there (see among others, Whitmore et al. 1993; 
Schweizer et al. 1996; Zepf et al. 1999; Alonso-Herrero 
et al. 2000, and 2001a - AAHOO, AAHOl; Scoville et al. 
2000 and references therein). This population of SSCs is 
not only inherent to LIRGs and ULIRGs, but it has also 
been identified in groups of galaxies (e.g., Gallagher et 
al. 2001) or even isolated barred or ringed galaxies (e.g., 
Alonso-Herrero, Ryder, & Knapen 2001b; Maoz et al. 2001 
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Table 1 
Sample and observations. 



Galaxy 


log LiR 


Dist 


Morphology 


Filters 


Linear 


Program 




(Le) 


(Mpc) 






Scale 




NGC 6808 


10.94 


46 


Isolated 


NIC3 F160W, F190N 


45 pc pixel" ^ 


7919 


NGC 5653 


11.01 


47 


Isolated 


NIC3 F160W, F190N 


45 pc pixel" ^ 


7919 


Zw 049.057 


11.22 


52 


Isolated, edge-on 


NIC2 F160W, F187N, F190N 


19 pc pixel" 


7219 


NGC 3256 


11.48 


37 


Advanced merger 


NIC2 F160W, F187N, F190N 


14 pc pixel" ^ 


7251 


NGC 1614 


11.62 


64 


Advanced merger 


NIC2 F160W, F187N, F190N 


24 pc pixel" ^ 


7218 


VV 114 


11.62 


80 


Interacting pair 


NIC2 F160W, NIC3 F187N, F190N 


29, 78 pc pixel" ^ 


7219 


NGC 6240 


11.82 


97 


Advanced merger 


NIC2 F160W, F187N, F190N 


36 pc pixel" ^ 


7219 


Arp 299* 


11.91 


42 


Interacting pair 


NIC2 F160W, F187N, F190N 


15 pc pixel" 


7218 



Note. — The infrared luminosity corresponds to the spectral range between 8 and 1000 /um. 
*The two components of Arp 299 are usually referred to as IC 694 (the Eastern component) and NGC 3690 (the Western 
component). 



and references therein). For a recent review on the sub- 
ject of SSCs in external galaxies, we refer the reader to 
Whitmorc (2000). 

While a lot of effort is being devoted to understanding 
the properties of SSCs, HST/NICMOS imaging has only 
recently revealed a population of bright H ii regions in two 
LIRGs, Arp 299 and NGC 1614 (AAHOO; AAHOl). A large 
fraction of these H II regions show luminosities in excess of 
that of 30 Doradus, the prototypical giant H il region. One 
of the main difficulties in quantifying the age of the stellar 
populations in LIRGs and interacting galaxies is break- 
ing the age-extinction degeneracy. This usually translates 
into only rough age estimates for SSCs (from photometric 
data) ranging between 5 and 900 Myr (see the recent re- 
view by Whitmore 2000). Hll regions, on the other hand, 
will highlight the youngest regions of star formation, with 
ages of < 5 — 10 Myr, as these are the lifetimes of the O and 
B stars required to ionize the gas. Clearly, understanding 
the properties of Hll regions and SSCs, and their relation 
at high spatial resolution will provide further insight into 
the nature of the star formation processes in LIRGs. 

In this paper we present a study of the detailed (tens 
to a few hundred parsecs) properties of the star forming 
regions (Hll regions and star clusters) of a sample of 8 
LIRGs. This paper is organized as follows. Section (2) 
describes the observations, data reduction and the pro- 
duction of the Hli region and star cluster catalogs. In 
Section (3) we establish the overall morphology of star 
forming regions in LIRGs and its relation with the dy- 
namical stage of the galaxy. In Section (4) the statistical 
properties of H ii regions in LIRGs are analyzed and com- 
pared with those of normal galaxies observed at compa- 
rable spatial resolutions. The spatial distribution of Hii 
regions and star clusters, their relative numbers and the 
age sequence are analyzed in Section (5). Our conclusions 
are presented in Section (6). 

iNote that although NGC 6808 is not technically a LIRG, its inf 
is included in our sample 



2. OBSERVATIONS 

The large amount of extinction routinely present in 
LIRGs, and in particular, the fact that active star form- 
ing regions are expected to contain non-negligible amounts 
of dust, prompted us to search the HST archive for in- 
frared observations. The obvious choice was narrow-band 
Paa and broad-band continuum imaging to identify Hll 
regions and star clusters respectively. This resulted in 
a sample of eight LIRGs. The sample covers a range of 
infrared luminosities - between logLiR — 10.94 Lq and 
logLiR = 11.82 Lq - as well as a variety of dynamical 
stages: isolated galaxies, close pairs of interacting galaxies 
and advanced mergers. The sample is presented in Ta- 
ble 1^ in increasing order of infrared luminosity. 

2.1. Data reduction 

The observations of the LIRGS analyzed in this paper 
were obtained as part of a variety of HST/mCMOS GTO 
and GO programs, listed in Table 1. The Paa-|-continuum 
images were taken with the NIC2 and NIC3 cameras (pixel 
size 0.076"pixel~"^ and 0.2" pixel^^, respectively) using the 
narrow-band filter F190N filter. At the distances of the 
LIRG sample, this filter (AA/A ^ 1%) contains the Paa 
emission line and the adjacent continuum at 1.90 /im. For 
the continuum subtraction, images through the F187N fil- 
ter were used, except for NGC 5653 and NGC 6808 (from 
Boker et al. 1999 survey) for which broad-band F160W fil- 
ter observations were employed instead. The field of view 
of the images is 19.5" x 19.5" and 51.2" x 51.2" for the 
NIC2 and NIC3 observations, respectively. All the contin- 
uum images used to identify the star clusters were observed 
through the F160W filter which represents a good com- 
promise between the better spatial resolution at shorter 
wavelengths and lower extinction at longer wavelengths. 

In Table 1 for each galaxy we list the infrared luminosity, 
distance (assuming Hq = 75kms^^ Mpc"^), large scale 
morphology, cameras and filters used, the corresponding 

luminosity (logLiR = 10.94 Lq) is close enough to the limit that it 
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NGC6240 X^^^=1.922^m, a=0.004/j,m 
W114 X^j^^=1.913/j,m, o-=0.003^im 




1.94 



A (/i.m) 



Fig. 1. — Observed Paa lines for VV 114 (dashed line) and NGC 6240 (solid line) depicted as Gaussian functions with an arbitrary flux, and 
the observed values of the wavelength and line width. The same emission lines are also displayed as if their observed wavelengths corresponded 
to that of the center of the filter. We also show the transmission curve of the F190N filter as the solid thick line. 



linear scale in parsec per pixel and the HST program num- 
ber of the observations. 

The images were reduced using nicred (McLeod et al. 
1997) routines, which involve subtraction of the first read- 
out, dark current subtraction on a readout-by-readout ba- 
sis, correction for linearity and cosmic ray rejection (using 
FULLFIT), and flat fielding. For the data reduced we used 
on-orbit dark and flatfield images. The angular resolution 
(FWHM) of the NIC2 and NIC3 images is approximately 
0.14" and 0.30", respectively, as measured from the point 
spread function (PSF) of stars in the images. 

2.2. Flux calibration and correction for off-center 
emission lines 

The flux calibration of the F160W, F187N and F190N 
images was performed using conversion factors based upon 
measurements of the standard star P330-E, taken during 
the Servicing Mission Observatory Verification (SMOV) 
program (M. J. Rieke, private communication 1999). 

The observed wavelength of Paa for the two most dis- 
tant galaxies in our sample, VV 114 and NGC 6240, is 
1.9131 jim and 1.9217 //m respectively. Figure 1 shows the 
observed Paa emission lines as Gaussian functions with 
line widths of cr 4.23 ± 0.66 x lO-^m for NGC 6240 
(Goldader et al. 1997) and cr = 3.2 ± 0.3 x IQ''^ nia for 
VV 114 (assuming that the line is unresolved and a 10% 
error, Doyon et al. 1995, and Goldader et al. 1997). For 
comparison we also display the same emission lines as if 
they had been observed at the center of the F190N fil- 
ter. As can be seen from this figure for both galaxies, 
the observed Paa emission lines are significantly displaced 
from the center of the narrow-band F190N filter. Assum- 
ing Gaussian line profiles, we estimate that the NICMOS 
filter may intercept only ~ 20% and ~ 2% of the full Paa 



line fiux for these galaxies, respectively. Given the large 
and uncertain correction implied to the true line flux, our 
results for these two galaxies are only qualitative. 

Throughout the paper we have converted the Paa lu- 
minosity into the more commonly used Ha luminosity, as- 
suming case B recombination (#^ — 8.7). 

^ ^ Paa 

2.3. Production of the Hll region and star cluster 
catalogs 

The Hii region catalogs were produced using the soft- 
ware REGION, kindly provided by Dr. C. H. Heller (see 
Pleuss, Heller, & Fricke 2000, and references therein for 
a detailed description). REGION is a semi-automated 
method to locate and compute statistics of Hii regions 
in an image, based on contouring, taking into account the 
local background. The lower limit for the size of an Hii 
region is set to 9 contiguous pixels, which for the distances 
of the galaxies in our sample corresponds to minimum lin- 
ear sizes (diameters) of between 42 and 234 pc. Each pixel 
must have an intensity above the local background of at 
least three times the rms noise of that local background 
(see Rand 1992 for more details on the criteria employed). 
After identifying the Hii regions, the program measures 
their position, size (area) and luminosity by subtracting 
the closest local background from the observed flux. For 
more details on the method and a discussion on the spa- 
tial resolution effects see Pleuss et al. (2000), and Alonso- 
Herrero & Knapen (2001). 

In Section 5 we will compare the number of detected 
Hii regions with the number of star clusters. To this 
end we need to determine the magnitude to which the 
H-h&nd luminosity distributions of the star clusters are 
complete. The positions of the star clusters were identified 
from the 1.6 /xm continuum images using the DAOFIND 
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task in iraf. As discussed in AAHOO for Arp 299, wo 
find that the star chisters in most of the systems studied 
here are located on a local background with a substantial 
gradient, and a large degree of crowding. Thus an auto- 
mated method to perform photometry as that provided 
by DAOFIND may yield to systematic errors, especially 
for the faintest star clusters. For this reason we have per- 
formed aperture photometry on the clusters of the three 
galaxies with the best spatial resolution (NGC 3256 and 
the two components of Arp 299) for which the results will 
be discussed in Section 5. 

The photometry was obtained through a relatively small 
aperture (diameter = 0.456"), so an aperture correction 
was applied to account for the total flux. The absolute H- 
band magnitudes of the star clusters in these two systems 
range up to Mh = —17 mag. The detection lower limit de- 
pends on the underlying galaxy background and the degree 
of crowding. We find that the luminosity distribution of 
the star clusters is complete down to Mh = —14 mag and 
Mh ^ -13.5 mag for NGC 3256 and Arp 299, respectively. 
This includes approximately 77 ± 5% of the detected clus- 
ters in both galaxies. The error accounts for uncertainties 
in the local background subtraction. For the luminosities 
and properties of the brightest star clusters in NGC 1614 
we refer the reader to AAHOl. Scoville et al. (2000) pre- 
sented photometry for the brightest near-infrared star clus- 
ters in Zw049.057, NGC 6240 and VV 114. As is the case 
for NGC 3256 and Arp 299, the cluster iJ-band magni- 
tudes are foimd to range up to Mh — —17.2 mag. 

3. MORPHOLOGY OF THE PAa LINE EMISSION IN LIRGS 

The top panels of Figure 2 arc the contours of the 
continuum-subtracted Paa line emission for the galaxies 
in our sample. We show the locations of the near-infrared 
star clusters and the emission peaks of the H ii regions as 
stars and open circles respectively (see Section 5). The 
bottom left panels of Figure 2 are the continuum emission 
images at 1.6 /im, whereas the bottom right panels are the 
continuum subtracted Paa images for the same fields of 
view. The images are displayed with the original HST ori- 
entation. These figures illustrate the variety of Paa line 
emission morphologies found in LIRGs. 

The isolated galaxies in our sample (NGC 6808, 
NGC 5653 and Zw 049.056) exhibit Hii regions along the 
spiral arms. NGC 6808, the LIRG in our sample with the 
lowest infrared luminosity shows the typical H ii region 
distribution of a spiral galaxy. Zw 049.056 appears as a 
highly inclined galaxy, and only shows faint Paa emission 
along the disk of the galaxy, perhaps as a result of the high 
extinction there (Scoville et al. 2000). NGC 5653 shows 
asymmetric Paa emission, with most and the brightest Hll 
regions located in one of the spiral arms. This galaxy has 
been classified as a lopsided galaxy (Rudnick, Rix, & Ken- 
nicutt 2000) . Lopsidedness is usually assumed to be a dy- 
namical indicator of a weak tidal interaction, that is, when 
the interaction does not destroy the disk of the galaxy. 
No bright nuclear emission is detected in NGC 5653 and 
NGC 6808. 

The Paa line emission of interacting/merging systems 
shows a more complex morphology, differing significantly 
from the stellar light distribution (continuum images). A 
common denominator for all these systems is the bright 
emission from the nuclei of the galaxies. The Paa nu- 



clear emission is usually compact, extending over scales of 
500 — 650 pc, except in NGC 1614 where the; Paa emission 
in the form of a ring of giant H ii regions with an approx- 
imate diameter of 600 pc (see AAHOl for more details). 

In the interacting/merging systems, besides the bright 
nuclear Paa emission, the recent star formation is also 
widely spread throughout the galaxies. In NGC 3256 and 
NGC 1614 the Paa emission extends along the spiral arms 
over at least 4kpc and 6kpc (that is, the field of view of 
the HSTfNlCMOS images for these galaxies). In the in- 
teracting systems Arp 299 and VV 114, Hii regions are 
detected throughout the galaxies, and at the interface re- 
gion. This is in good agreement with findings for ULIRGs 
in Murphy et al. (2001). The Paa morphology in the 
merging/interacting systems clearly shows that the effects 
of extreme star formation can propagate throughout the 
galaxies and not only in the nuclei. Recently Bekki & 
Couch (2001) have run stellar and gas numerical simula- 
tions of galaxy mergers to find the most promising sites 
for formation of SSCs. They conclude that these SSCs are 
not only formed in the central regions, but also between 
the two interacting galaxies. 

NGC 6240 only shows Paa emission originating in both 
nuclei (separated by only 0.8 kpc). This is consistent 
with recent HST/WFPC2 Ha imaging (Gerssen et al. 
2001). The larger scale filamentary diS^use morphology 
in NGC 6240 has not been detected by the NICMOS ob- 
servations, probably as a result of the Paa line falling at 
the edge of the narrow-band filter (see Section 2.2). 

As most of the infrared luminosity in LIRGs is believed 
to be due to dust heated in regions of strong star for- 
mation, the mid-infrared emission in LIRGs would be ex- 
pected to trace the Paa emission. Recently Soifer et al. 
(2001) have obtained high resolution mid-infrared obser- 
vations of LIRGs and reported that a substantial fraction 
(> 50%) of the infrared luminosity is generated in re- 
gions of sizes ranging from 100 pc to a 1 kpc, with the 
nuclear starbursts often dominating the mid-infrared lu- 
minosity. They also concluded that the 12 /xm emission 
traces the current sites of dusty luminous star formation. 
This can be confirmed by the good correspondence be- 
tween the mid-infrared and the Paa emission morpholo- 
gies of the galaxies in common with our study: Arp 299, 
VV 114 and NGC 1614. 

In Table 2 we list the total Ha luminosity from the de- 
tected H II regions, along with the percentage of the Paa 
emission originating in the nuclei of galaxies. These frac- 
tions should be considered only as lower limits as the com- 
pact nuclear regions are usually more reddened than the 
regions where the more diffuse H ii regions are located (see 
e.g., AAHOO; Lfpari et al. 2000; Scoville et al. 2000; see 
also Murphy et al. 2001 for ULIRGs). Nevertheless it 
is apparent that the nuclear to total emission ratios vary 
significantly from system to system, and show similarities 
with the same ratios computed for the infrared luminosi- 
ties (Soifer et al. 2001). In the isolated systems most of 
the star formation is occurring in the disks of the galaxies, 
whereas in the interacting/merging systems the location of 
the star forming regions depends on factors such as the age 
of the interaction process and the physical properties of the 
interacting galaxies (among others, the initial gas content 
and mass ratios of the interacting galaxies, see e.g., Bekki 
& Couch 2001). Table 2 also gives the ratios between the 
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Table 2 

Line emission and Hii region properties in LIRGs. 





FOV nf Parr 






\o(>:( Ltd /LCHrf")* A 






(erg s"^) 






NGC 6808 


11 kpc X 11 kpc 


41.38 


~0 


3.1 


NGC 5653 


11 kpc X 11 kpc 


41.72 


~0 


2.9 


Zw 049.057 


4.9 kpc X 4.9 kpc 


41.21 




3.6 


NGC 3256 


3.5 kpc X 3.5 kpc 


42.11 


0.41 


3.0 


NGC 1614 


6.4 kpc X 6.4 kpc 


42.60 


0.65 


2.5 


VV 114 


9.3 kpc X 9.3 kpc 


42.55* 


0.37 


2.7* 


IC 694 


3.8 kpc X 3.8 kpc 


41.95 


0.58 


2.6 


NGC 3690 


3.8 kpc X 3.8 kpc 


42.16 


0.18 


3.1 


NGC 6240 


9.2 kpc X 9.2 kpc 


43.19* 


~ 1* 


2.2* 



Note. "FOV of Paa" is the area imaged through this emission hnc. 
* Because of the large correction needed to account for the Paa emission hne flux 
missing in the F190N filter (Section 2.2), these values are uncertain. 



Fig. 2a. — Bottom left panel: Grey-scale continuum image at 1.6 fim (on a logarithmic scale) of NGC 6808. Bottom right panel: Grey-scale 
continuum subtracted Paa image (on a logarithmic scale). Upper panel: The contours (on a linear scale) axe the continuum-subtracted Paa 
emission. The stars (blue stars in the on-line color version of Figure 2) are the positions of the identified near-infrared star clusters, whereas 
the open circles (red circles in the on-line color version of this figure) are the emission peaks of the identified H il regions. For reference, we 
mark the position of the nucleus with an 'N'. 



Fig. 2b. — Continued: Same as Figure 2a but for NGC 5653. 



Fig. 2c. — Continued: The left panel is the continuum image of Zw 049.057. The right panel are the contours of the Paa emission, with 
only the positions of star clusters shown. 



Table 3 

Statistical properties of Hii regions. 



Galaxy 


logL(Ha)br 


log < L(Ha) >3 


logL(Ha)mcdian 


SlZGlarge 


< size >3 


SlZ6median 




(erg s"^) 


(erg s"^) 


(erg s-i) 


(pc) 


(pc) 


(pc) 


LIRGs 


NGC 3256 


41.53 


41.24 


39.51 


308 


253 


80 


IC 694 


41.71 


41.29 


39.56 


448 


298 


82 


NGC 3690 


41.43 


41.33 


39.72 


346 


320 


94 


NGC 1614 


41.01 


41.00 


39.89" 


210" 


205" 


148" 


NGC 6808 


40.18 


40.11 


39.22 


461 


423 


201 


NGC 5653 


41.02 


40.74 


39.49 


551 


468 


199 


VV 114 


41.94** 


41.79*' 


40.52*" 


920 


842 


382 


Comparison Normal Galaxies 


NGC 4192 


39.52 


39.39 


38.52 


158 


147 


76 


NGC 4389 


39.44 


39.36 


38.70 


151 


134 


74 


IC 750 


39.89 


39.73 


38.62 


191 


184 


72 


NGC 4102 


40.48 


40.32 


38.32 


263 


223 


80 



Note. — " Statistics excluded the bright H ii regions in the ring of star formation, as blending effects are present. 
' Uncertain because of the large correction needed to account for the Paa emission line flux missing in the F190N 
filter (Section 2.2). 




NGC3256-H NGC3256-Paa 




5 -5 5 -5 



arcsec arcsec 

Fig. 2d. — Continued: Same as Figure 2a but for NGC 3256. 
Fig. 2e. — Continued: Same as Figure 2a but for NGC 1614. 

Fig. 2f. — Continued: Same as Figure 2a but for VV 114. The NIC2 F160W continuum image was degraded to the resolution of the NIC3 
F190N image before the identification of the star clusters (as shown here). The nuclei of the system are marked with the symbols 'NE', 'NW' 
and 'W (following Soifer et al. 2001 notation). 
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Fig. 2g. — Continued: Same as Figure 2a but for IC 694, the Eastern component of Arp 299. 



Fig. 2h. — Continued: Same as Figure 2a but for NGC 3690, the Western component of Arp 299. 



Fig. 2i. — Continued: The left panel in the continuum image of NGC 6240. The right panel are the contours of the Paa emission with the 
positions of the clusters and H ii regions marked as in Figure 2a. 



infrared luminosity and the total Ha luminosity for the 
LIRGs in our sample. These ratios are found to be sim- 
ilar for all the galaxies except for Zw 049.057 where the 
Paa emission is probably affected by high extinction, and 
NGC 6240 for which the Paa luminosities are uncertain as 
discussed in Section 2.2. 

4. STATISTICAL PROPERTIES OF GIANT H II REGIONS IN 

LIRGS 

4.1. Luminosities, Sizes and Luminosity Functions 

We have measured the luminosity and diameter of the 
identified H ii regions in our sample of LIRGs. For each 

galaxy wc determined the Ha luminosity for the brightest 
H II region, for the three brightest H ii regions (first-ranked 
H II regions) and the median H II region, as well as the di- 
ameter for the largest Hii region (which is not necessar- 
ily the brightest Hii region), for the first-ranked Hll re- 
gions and the median H ii regions. The results are summa- 
rized in the first part of Table 3, where galaxies are sorted 
in order of decreasing spatial resolution. As discussed in 
Alonso-Herrero & Knapen (2001) and others, even at HST 
spatial resolutions the problem of defining an Hll region 
remains. For instance, at increasing higher spatial resolu- 
tions different emission peaks within an H il region may be 
identified as individual Hii regions thus cataloging them 
as smaller and less luminous Hii regions. Conversely, as 
the spatial resolution becomes poorer, blending effects will 
become more frequent. This effect is clearly seen in the in- 
creasing median sizes of H II regions in the LIRG sample 
for decreasing spatial resolutions (Table 3). 

All the LIRGs in our sample contain a significant num- 
ber of exceptionally bright Hii regions (i.e., the median 
values of the H II region luminosity distribution) with Ha 
luminosities comparable to that of the giant H ii region 30 
Doradus (log L(Ha) ~ 39.70ergs"i). The first-ranked H ii 
regions in the interacting/merging galaxies usually coin- 
cident with the nuclei of the galaxies - are between one and 
two orders of magnitude brighter than the median H II re- 
gions. As can be seen from Table 3 the luminosities of the 
brightest and first-ranked H II of the two isolated galax- 
ies, even at poorer spatial resolutions, are markedly less 
luminous than their counterparts in interacting/merging 
systems. Unfortunately the small number of LIRGs in 
our sample does not allow us to establish whether isolated 
LIRGs tend to have fainter H ii regions than interacting 
LIRGs. 

The statistical properties of H ii regions are represented 
with LFs, usually fitted by a power law: 

dN{L)=AL"dL (1) 

where a is the index of the power law with values of 
a = —2.0 ± 0.5 for extragalactic Hii regions (Kcnnicutt, 
Edgar, & Hodge 1989). We have constructed Hii re- 
gion LFs for the two LIRGs in our sample with a suffi- 



ciently high number of detected Hll regions: Arp 299 and 
NGC 3256 (shown in Figure 3). We have fitted the slopes 
of the LFs excluding the nuclei and the regions of strong 
star formation activity in Arp 299 (as they may be blends 
of a few regions) and found values of a = —1.77 ± 0.10 
and a = -1.81 ± 0.16 for NGC 3256 and Arp 299, respec- 
tively. If we include all Hii regions in the fit, the values 
of the slopes are a = — 1.68 ± 0.08 and a = — 1.66 ± 0.05 
respectively. These slopes are within the values fitted for 
H II regions in the disks and circumnuclear regions of nor- 
mal galaxies (e.g., Kennicutt et al. 1989), and those found 
from fitting the luminosity functions of star clusters (e.g., 
Zepf et al. 1999; Whitmore 2000). 

4.2. Comparison with normal galaxies 

A number of studies have analyzed the statistical prop- 
erties of Hii regions (luminosities, sizes, density numbers) 
in the disks of normal spiral galaxies in terms of the mor- 
phological type of the galaxy and the arm-interarm re- 
gions (see for instance the classical study of Kennicutt et 
al. 1989). Because of the effect of the spatial resolution on 
the properties of H ii regions, we chose to compare the sta- 
tistical properties of the H II regions in LIRGs with those of 
the circumnuclear H II regions of normal galaxies (as de- 
rived from _ff5T/NICM0S Paa observations) studied by 
Alonso-Herrero & Knapen (2001). From this work, we 
have selected galaxies at the distance of the Virgo Cluster 
(Dist ~ 17Mpc), for which the NIC3 observations provide 
a resolution of 16pc pixel" ^, similar to that of NGC 3256, 
NGC 3690 and IC 694. 

In the second part of Table 3 we summarize the statis- 
tical properties of H II regions of the normal galaxies with 
the highest degree of central star forming activity (as mea- 
sured from the central kpc Ha luminosity): NGC 4192, 
NGC 4389, IC 750, and NGC 4102. 

Figure 4 shows histograms comparing the Ha luminosi- 
ties of Hii regions in LIRGs and normal galaxies. We 
have divided the normal galaxies into two groups, isolated 
galaxies (NGC 4192, NGC 4389 and NGC 4102) and the 
interacting galaxy IC 750. In this figure we also indicate 
the Ha luminosity of 30 Doradus, the prototypical giant 
H II region. These histograms not only show that giant H II 
regions are more common in LIRGs than in normal galax- 
ies - even when comparing with the interacting system 
IC 750 - but also that the median Hii regions in LIRGs 
are at least an order of magnitude brighter than in normal 
galaxies. 

The large population of giant Hii regions in LIRGs - 
with a significant fraction being more luminous than 30 
Doradus - is unprecedented in normal galaxies. There are 
a number of possible explanations. We may just be seeing 
the extended tail of the H ii region LF which translates into 
a larger number of the more massive ionizing star clusters, 
a extremely young population of star forming regions in 
LIRGs, or perhaps the bright Hii regions in LIRGs may 
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Fig. 3. — Hll region LFs for NGC 3256 and Arp 299. The solid lines are the power law fits excluding the nuclei of the galaxies, and the 
regions with strong star formation activity in Arp 299 (C and C). For comparison the dashed lines represent the fits including all the sources. 



represent aggregations of normal Hll regions. The typical 
(median) diameters for H ii regions in normal galaxies and 
LIRGs observed at the same spatial resolution are compa- 
rable (see Table 3). This seems to argue against the last 
possibility. Even though the first ranked Hll regions in 
LIRGs are typically a factor of two larger than in normal 
galaxies, the luminosities of the first-ranked (as well as 
median) H li regions in LIRGs would require aggregations 
of ten "normal" Hll regions. Other possibilities are then 
evolutionary effects and the upper limit to the mass of the 
ionizing clusters. 

As can be seen from the simulations of H ii region LFs 
by Gey & Clark (1998) the age of the stellar population is 
one of the factors determining the high luminosity end of 
the H II region LF. For a population of H ii regions formed 
in a single burst, the high luminosity end of the LF as well 
as the slope vary as the population ages. In the case of 
continuous creation of H II regions the high luminosity end 
of the LF remains approximately constant. The youth of 
a population of coeval Hii regions could be a possibility 
to explain the excess of bright Hii regions in LIRGs. As 
we shall see in the next section however, the Hii regions 
in LIRGs show a range of ages, which means that the ex- 
treme luminosities of the H II regions in LIRGs are not due 
to the fact that they are all very young. 

Since the fitted slopes of the Hii region LFs of two 
LIRGs are consistent with the values found in normal 
galaxies, the extraordinary luminosities of H II regions in 
LIRGs are simply reflecting a larger number of massive 
star clusters in LIRGs, rather than an anomalous stellar 
mass distribution. As discussed in Bekki & Couch (2001), 
the physical conditions (high gas pressure and density of 
the interstellar medium) necessary for the formation of 
super star clusters can only be achieved in systems under- 
going a rapid transfer of gas to the central regions, that 
is, interacting galaxies. It is unlikely that these conditions 
are met in isolated disk galaxies. 

5. GIANT HII REGIONS AND THEIR RELATION TO SSCS 



5.1. Spatial distribution of the Hii regions and 
near-infrared star clusters 

Figure 2 shows that the emission peaks of the Hii re- 
gions are generally not spatially coincident with the loca- 
tion of the near- infrared star clusters. To quantify this, 
we have cross-correlated the positions of the H ii regions 
and the near- infrared clusters. We counted a coincidence 
when the spatial separation between an H ll region and a 
cluster was equal or less than 0.15"-0.11" (for the NIC2 
images) and 0.40"-0.30 "(for the NIC3 images), that is, a 
separation of 2 and 1.5 pixels respectively. 

The results from the cross-correlation are presented in 
the first part of Table 4 where we list the number of de- 
tected Hii regions, and star clusters, together with the 
number of coincidences. As was clear from Figure 2, 
the number of coincidences is relatively small. The in- 
creasing number of coincidences (intermediate age H ll re- 
gions/clusters) for NGC 5653, NGC 6808 and VV 114 may 
be the result of the poorer spatial resolutions of the NIC3 
images, which implies that a coincidence is counted when 
the spatial separations are 68 — 117 pc (90 — 156 pc for the 
2 pixel separation). For comparison, the spatial separa- 
tions for the galaxies observed with the NIC2 camera are 
21 — 38pc (or 28 — 50 pc if we use the 2 pixel separation). 
This lack of spatial coincidence between Hii regions and 
clusters has also been found in barred and ringed galaxies 
with strong star formation activity (Alonso-Herrero et al. 
2001b; Maoz et al. 2001). 

Theoretical arguments may offer an explanation for the 
apparent offsets between the location of the SSCs and the 
H II regions in LIRGs. For instance. Tan & McKce (2000) 
have proposed a model in which young clusters within 
clouds with masses 5 x 10^ — 5 x 10"'' M0 can survive for 
up to 3Myr and allow for the formation of gravitation- 
ally bound clusters (which in turn may evolve into glob- 
ular clusters) in the presence of vigorous feedback (Wolf 
Rayet winds and supernovae). According to their model, 
after approximately 3Myr, the gas will be disrupted and 
the star cluster will emerge from its cocoon, and will be 
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Fig. 4. — Histograms comparing the Ha luminosities of the Hll regions detected in two LIRGs (Arp 299 and NGC 3256) in the bottom 
panel, normal galaxies (NGC 4192, NGC 4389 and NGC 4102) in the middle panel and the interacting galaxy IC 750, in the top panel. All 
these galaxies have been observed with comparable spatial resolutions. The arrow indicates the Ha luminosity of the giant H ii region 30 
Doradus. 
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visible. The youngest SSCs are still embedded in the 
(probably dusty) H II regions, and demonstrates that gi- 
ant H II regions offer us an unprecedented view to the truly 
youngest (high-mass) star formation activity. The SSCs on 
the other hand show a broader range of ages as we shall 
show in the next section. 

There are also Galactic examples of this situation in ob- 
scured giant Hll regions, where a variety of complex mor- 
phologies including nebular emission, regions of high and 
variable extinction and centrally concentrated star clusters 
are seen (see Blum et al. 2001 for a census, and references 
therein). In particular these authors have established an 
evolutionary sequence for their sample of giant H Ii regions 
(see their figure 1), based on the morphology and spectro- 
scopic ages of the stars. In the youngest H ii regions the 
cluster stars are still veiled by the hot dust from their birth 
cocoons; as the H ii region ages, the clusters become more 
revealed as the gas will be dispersed by energetic winds and 
radiation pressure from the hot stars, until finally most of 
the gas will have dispersed. 

The lack of spatial correspondence between the near- 
infrared SSCs and the bright H II regions in LIRGs is sug- 
gestive of a similar evolutionary sequence as that seen in 
obscured Galactic giant Hii regions. To first order, ig- 
noring possible effects of extinction and the observational 
detection limits, those Hll regions with no near-infrared 
cluster counterpart will tend to be the youngest exam- 
ples of star forming regions in which the gas is being 
ionized by the most massive stars, while the star clus- 
ters will still be embedded; we denote them as "young 
H II regions/clusters" . They account for between 20% and 
40% of all the detected sources (star clusters and H ii re- 
gions) in LIRGs (see Table 4). The coincidences corre- 
spond to the case of an evolved Hii region in which the 
near-infrared cluster is becoming "visible", that is, the 
time period in which the cluster has begun to emerge from 
the natal dust. At the same time the first red supergiants 
will start appearing so the cluster will become bright in the 
near-infrared. The "intermediate Hii regions/clusters" ac- 
count for up to 30% of the cluster -|- H ii region population 
in LIRGs. We also find that the number of coincidences 
tends to increase for the more luminous Hii regions. For 
instance, in NGC 3256 we find that ~ 55% of the high 
luminosity Hll regions have a new-infrared cluster coun- 
terpart (see next section). And finally, the near- infrared 
star clusters with no H ii region counterpart represent the 
case where the most massive (ionizing) stars have disap- 
peared so there will be no detectable gas emission; these 
will be referred to as "old clusters" . We estimate that they 
contribute between 56% and 72% of the detected cluster 
population of LIRGs. 

5.2. Modelling of the Hii region to near-infrared star 

cluster age sequence 

The age/coincidence argument can be developed fur- 
ther by using outputs of evolutionary synthesis models 
and the relative numbers of Hii regions, coincidences and 
near-infrared star clusters. Obviously for a given galaxy, 
the number of identified sources depends on the detection 
threshold of the current observations. It will be easier 
to present this argument using the three galaxies with the 
best spatial resolutions: NGC 3256, IC 694 and NGC 3690. 
The argument however will be valid for any galaxy by 



changing the detection limits. 

In Section 2.3 we established that the i7-band lumi- 
nosity functions of the star clusters in NGC 3256 and 
Arp 299 are complete down to absolute magnitudes of 
Mh = —14 mag and Mh = —13.5 mag, respectively. The 
H II region luminosity functions are approximately com- 
plete down to logL(HQ:) = 38.8 — 39 erg s~"^ (previous sec- 
tion), which is equivalent to a number of ionizing photons 
of logA^Ly = 50.7 — 50.9s~^. We can use these detec- 
tion limits to compute the fractions of young, intermedi- 
ate and old populations for the complete distributions of 
NGC 3256, IC 694 and NGC 3690 (second part of Ta- 
ble 4). Using the complete distributions increases the Hii 
region fraction by ~ 5%, whereas the old cluster fraction 
drops by ~ 5%, and the coincidence fraction remains ap- 
proximately constant. We will make use of evolutionary 
synthesis models to establish an age sequence, as well as 
to reproduce the fractions of young, intermediate and old 
clusters in these systems. Our choice for the evolution- 
ary synthesis models will be Starburst99 (Leitherer et al. 
1999) and those of Rieke et al. (1993). 

In Figure 5 we show outputs of Starburst99 for two dif- 
ferent cluster masses: 10^ Mq and 10^ Mq, a Salpeter IMF 
(with lower and upper mass cutoffs of 1 M0 and 100 M0 re- 
spectively), instantaneous star formation and solar metal- 
licity. The solid lines represent the time evolution of the 
number of ionizing photons (thin solid line) and the H- 
band absolute magnitude (thick solid line). For simplicity 
we only show the detection thresholds for NGC 3256 as 
dashed-dotted lines at Mh = —14 mag (scale on the right 
hand side of the diagram) and logiVLy = 50.7 s~^ (scale 
on the left hand side of the diagram) . For a 10^ Mq clus- 
ter (lower panel of Figure 5), and for ages up to 2± 7Myr 
we will observe a coincidence between an Hll region and 
a near-infrared star cluster (i.e., when both A^Ly and Mh 
are above their respectively detection thresholds, shown 
as the closely hatched area), assuming no extinction (see 
Section 5.4.1). After 7Myr and for at least lOOMyr, we 
will be able to detect star clusters (grey area), but no co- 
incidences, since A^Ly has fallen below its threshold. 

The case of a Starburst99 10^ Mq cluster is displayed 
in the upper panel of Figure 5. The detection limits in 
NGC 3256 will allow us to identify Hii regions for up to 
5 Myr (shown as a lightly hatched area) and near-infrared 
star clusters for ages 7—12.5 Myr (shown as the grey color 
area). After approximately 12 — 13 Myr, the star clusters 
will be too faint to be detected with the current detection 
threshold. Note that for clusters of 10^ Mq mass there 
will be no age for which we could detect both an Hii re- 
gion and a near- infrared star cluster (that is, there will be 
no coincidences). Within the present detection limits, for 
NGC 3256 and Arp 299 star clusters with masses of less 
than 10^ Mq will not be detected as either star clusters or 
H II regions, as at no time during their evolution will they 
be brighter than Mh = —14 mag or logA^Ly = 50.9 s~^, 
respectively. 

For a more extended period of star formation we have 
used outputs of Rieke et al. (1993) models with a Gaussian 

star formation rate of 5 Myr FWHM. The results are pre- 
sented in Figure 6 for two clusters with masses 5x10^ Mq 
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Fic. 5. The solid lines are outputs of Staiburst99 (Lcithcrcr et al. 1999) showing the time evolution of the number of ionizing photons 
(log A^Ly, thin solid line, scale on the left hand side) and absolute /f-band magnitude {Mh, thick solid line, scale on the right hand side) for a 
10® M0 cluster (bottom panel) and a 10^ M.q cluster (top panel), for an instantaneous burst, a Salpeter IMF (between 1 and 100 Mq) and solar 
metallicity. The dashed-dotted lines represent the approximate detection thresholds for NGC 3256: Mh = — 14mag and logA^Ly = 50.7s~^. 
An H 11 region with no cluster counterpart will be the case when log A^Ly is above the detection threshold, but Mj{ is not detectable yet 
(lightly hatched area) . There will be a coincidence between an H ii region and a near-infrared star cluster when both Nj^y and Mh are above 
their respectively detection thresholds (closely hatched area). We will only observe a star cluster when logATLy < 50.7 s~^, but Mh is still 
observable (grey area). 
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Fig. 6. — Same as Figure 5, but for outputs of Rieke et al. (1993) models, using a truncated Salpeter IMF (between 0.1 and 80 M©), a 
Gaussian burst of 5Myr duration and a cluster with 5 x 10^ Mq mass (lower panel) and 5 x 10'' Mq (upper panel). 
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and 5 x 10^ Mq for a truncated Salpeter IMF^. Figure 6 is 
coded in the same way as Figure 5. It is clear that the dif- 
ferent star formation rates will produce different fractions 
of Hii regions, coincidences and star cluster detections. 
For instance, in the case of a 5 x 10^ M0 cluster formed in 
a Gaussian burst, we will detect Hii regions alone between 
3 and 9Myr, before the 7?-band luminosity of the cluster 
becomes detectable. In the case of the 5 x 10^ Mq cluster, 
because the star formation period is more extended than 
in the instantaneous burst, it is still possible to have coin- 
cidences for a short period of time, between 9 and 12Myr. 

We now need to assume a mass distribution for the star 
clusters. Elmegreen & Efrcmov (1997) proposed that all 
types of clusters (young star clusters, globular clusters, 
open clusters and associations) form with nearly constant 
efficiency in gas clouds. This universal mechanism applies 
then to all regions regardless of epoch or geometry. Thus 
the mass distribution of the star clusters in NGC 3256 and 
Arp 299 can be represented as a power law: 

n{M) AM oc M"" AM (2) 

where a is the slope of the distribution. From star cluster 
and Hii region LFs, it is found that the slope has values 
of a = 1.7 — 2 (see Elmegreen & Efremov 1997 and ref- 
erences therein). In particular, Zepf et al. (1999) have 
found a slope of a = 1.8 for the optical LF of the young 
star clusters in NGC 3256. 

The next step is to assume an age distribution for the 
star clusters. From optical spectroscopy of the H II regions 
of Arp 299, we measured a range of equivalent widths of 
Ha (a similar situation is found in NGC 3256, see Li'pari et 
al. 2000) which implies that the star clusters/Hii regions 
are not coeval. We will assume that clusters are formed at 
a constant rate, so there is a range of ages. Using Equa- 
tion (2) for the mass distribution of the star clusters, and 
the evolutionary synthesis models presented above, we can 
obtain an estimate the fractions of young, intermediate 
and old clusters. For the two galaxies in consideration 
wc will only need to take into account star clusters with 
masses between ~ 5 x 10" Mq and 10'^ Mq. Less massive 
clusters will not be observed with the present detection 
threshold, whereas the upper limit of the iJ-band lumi- 
nosity distribution [Mh — —17.2 mag) indicates that there 
are no star clusters more massive than lO'' Mq. The latter 
seems to be true for all the interacting/ merging LIRGs in 
our sample. 

5.3. Comparison with observations 

The first result is that the relative numbers of young 

H II regions/clusters, intermediate and old clusters depend 
only slightly on the slope of the mass distribution (for the 
range a = 1.7 — 2). The main factors determining these 
fractions are the age distribution of the star clusters (that 
is the age of the oldest clusters), and the adopted form 
for the IMF and star formation rate. See Sections 5.4 and 
5.5 for a discussion on the effects of extinction and cluster 
destruction. 

In Table 5 we summarize the results for a cluster mass 
distribution with a slope of a = 1.7. The first part of the 

table arc the results for an instantaneous burst, whereas 
the second part are those for a more extended period 

'^<f>{m)Am oc mT^'^^Am between 1 and 80 Mq, and (/)(m)dm oc m~ 



of star formation. The model predictions can be com- 
pared with the fractions for the complete distributions of 
NGC 3256 and Arp 299 given in the second part of Table 4. 

Inspection of Figures 5 and 6, and taking into account 
that less massive clusters are more numerous than more 
massive clusters, shows that the more extended star for- 
mation activity of Rieke et al. (1993) models will always 
produce more coincidences than observed, regardless of the 
maximum age of the star clusters. This probably indicates 
that the instantaneous star formation is a better choice to 
account for the observed properties of star clusters. This 
seems to be the case for local star clusters where the age 
spread is at most just a few million years (e.g., Luhman 
et al. 1998; Palla & Stahler 2000). This is not true in 
general for the nuclei of interacting galaxies, and regions 
of high star formation activity (see AAHOO and AAHOl) 
where we have found that more extended periods of star 
formation are required to explain the observed properties. 

Figures 5 and 6 also illustrate the dependence of the 
cluster luminosity with the model assiuiiptions, which in 
turn translates into different mass estimates from the ob- 
served absolute i?-band magnitudes. Using the instanta- 
neous star formation and a Salpeter IMF between 1 and 
100 Mq we can constrain the masses of the identified star 
clusters in NGC 3256 and Arp 299 between ~ 5 x lO" M© 
and ~ 10^ Mq. The high end of our photometric masses 
is in good agreement with the dynamical masses measured 
by Mengcl ct al. (2002) for compact young star clusters in 
the Antennae system. 

Another interesting result from our simulations is that 
the observed relative fractions of young (H 11 regions) and 
old (near-infrared star clusters) populations provide some 
information on the age spread of the last epoch of star 
formation. From the statistical point of view our simula- 
tions show that most of the detected stellar population 
in Arp 299 and NGC 3256 shows an age spread of be- 
tween 20Myr and 40Myr, as otherwise we should have 
measurcid a higher fraction of old clusters (e.g., for an age 
distribution of up to 100 Myr this faction is ~ 80%) . Zepf 
et al. (1999) also deduced young ages for the NGC 3256 
star clusters detected in the optical based on the lack of 
a strong color-luminosity relation. The youth of the de- 
tected clusters does not exclude the presence of older clus- 
ters possibly created in this or previous episodes of star 
formation. These clusters however cannot be identified 
with the present detection threshold. 

5.4. Limitations 

In this section we discuss some limitations of the com- 
parison between the observed fractions of H 11 regions, old 
clusters and the model predictions. 

5.4.1. The effects of extinction 

So far we have ignored the effects of extinction in our 

comparison of observations and models. From detailed 
studies of obscured Galactic giant H 11 regions there is ev- 
idence that the extinction to the stars is highly variable 
with values between Ak = and Ak = 4.5 mag (Blum, 
DamineU, & Conti 1999; Blum, Conti, & DamineU 2000) 
for the very young H 11 regions. If high extinction is present 
in the Hii regions of LIRGs, then the predicted fractions 

^dm between 0.1 and 1 Mq. 
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Table 4 

Observed Numbers and Fractions of Hii regions and near-infrared star clusters. 



Galaxy 


No. Hii 


No. 


No. star 


Young H 11 


Intermediate 


Old 




regions 


Coincidences clusters 


regions 


clusters/H ii regions 


clusters 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC 3256 


129 


27-18 


221 


32 - 33% 


8-5% 


60 - 62% 


IC 694 


58 


10-7 


131 


27 - 28% 


6-4% 


67 - 68% 


NGC 3690 


75 


25 - 19 


210 


19 - 21% 


10 - 7% 


71 - 72% 


NGC 1614 


ring(20) + 30 


16 - 11 


99 


26 - 28% 


12 - 8% 


62 - 64% 


NGC 6808 


81 


18 - 11 


122 


34 - 36% 


10 - 6% 


56 - 58% 


NGC 5653 


71 


24 - 17 


144 


24 - 27% 


14 - 9% 


62 - 64% 


VV 114 


35 


20-15 


52 


22 - 28% 


30 - 21% 


48 - 51% 


Complete distributions 


NGC 3256 


128 


27-18 


170 


37 - 39% 


10 - 7% 


53 - 54% 


IC 694 


57 


9-6 


99 


33 - 34% 


6-4% 


61 - 62% 


NGC 3690 


73 


22 - 16 


158 


25 - 27% 


10 - 7% 


65 - 66% 



Note. — The range in numbers in columns (3), (5), (6) and (7) correspond to the criteria used for counting a 
coincidence between an Hii region and a cluster: a separation of 2 — 1.5 pixels 



Table 5 

Predicted fractions of Hii regions, coincidences and old clusters for NGC 3256 and 

Arp 299 (no destruction of clusters). 



Galaxy 


Age 


Young H II 


Intermediate 


Old 




Distribution 


regions / clusters 


clusters/H 11 regions 


clusters 


Instantaneous burst 


NGC 3256 


0-10 Myr 


50% 


14% 


36% 




0-20 Myr 


36% 


10% 


54% 




0-40 Myr 


28% 


8% 


64% 




- 100 Myr 


17% 


5% 


78% 


Arp 299 


0-10 Myr 


50% 


14% 


36% 




0-20 Myr 


31% 


9% 


60% 




0-40 Myr 


25% 


7% 


68% 




- 100 Myr 


16% 


4% 


79% 


5 Myr Gaussian burst 


NGC 3256 


0-10 Myr 


80% 


20% 


0% 


and Arp 299 


0-20 Myr 


34% 


23% 


43% 




0-40 Myr 


29 - 27% 


19 - 18% 


52 - 55% 




- 100 Myr 


20 - 19% 


13 - 12% 


67 - 69% 
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of Hii regions, coincidences and older star clusters will 
be affected. Although we have no way derive the extinc- 
tions from the current Paa observations, we can obtain 
some estimates from optical spectroscopy. AAHOO and 
Lipari et al. (2000) have measured extinctions of up to 
Ay — 4.2 mag using the Balmer decrement for a few Hii 
regions in Arp 299 and NGC 3256 respectively. This is 
equivalent to extinctions at the observed wavelength of 
Paa of up to Apaa 0.5 mag. The ages of these Hii re- 
gions from the observed equivalent widths of Ha are of 
between 3 and 6Myr using Leitherer et al. (1999) models. 

Although the measured extinctions are not high enough 
to compromise the detection of H II regions at near- infrared 
wavelengths, we cannot rule out the possibility that very 
young Hii regions (< 3Myr) suffer from elevated extinc- 
tions, even at near-infrared wavelengths. If this were the 
case, we would be missing the youngest H ii regions, and 
hence the observed fraction of H ii regions and coincidences 
will be lower limits. This would in turn translate in even 
younger age distribution of the detected population of star 
clusters. 

5.4.2. The effects of destruction of clusters 

The observed luminosity and mass functions of old glob- 
ular clusters around galaxies and young star clusters ob- 
served in interacting galaxies appear to be distinctively 
different, with the former having a log-normal shape, and 
the latter a power law form with no evidence for a turnover 
(see for instance Elmcgrcen & Efremov 1997; Whitmore 
et al. 1999; Zepf et al. 1999). This poses an interesting 
problem if the clusters observed in galaxies are younger 
versions of today's globular clusters. One of the proposed 
solutions is that the mass distribution of old globular clus- 
ters was initially a power law that was later modified by 
selective destruction of low mass star clusters to become a 
log-normal like mass function (see Fall & Zhang 2001 for 
a detailed discussion, and also Whitmore 2000). 

Although including a detailed treatment of the problem 
of cluster destruction is beyond the scope of this paper, 
we can attempt to see its effects on our calculations. If 
clusters in Arp 299 and NGC 3256 have been forming at 
a constant rate for approximately 100 Myr, then ~ 50% 
of all clusters with masses above 5 x 10'' should be 
destroyed during that period to account for the observed 
fraction of star clusters. Note that Zepf et al. (1999) 
proposed selective destruction of low mass clusters as one 
possibility to account for the observed optical colors and 
luminosities of clusters in NGC 3256 (the other possibility 
was a very young age for the clusters). 

5.4.3. Paa. shells 

Finally we consider the effects of the formation of Paa 
shells on the measured number of coincidences. In the An- 
tennae galaxy there is evidence that many of the slightly 
older, more massive clusters (i.e., 5-10 Myr) have blown 
large Ha shells around themselves, hence there is no longer 
a good correspondence between the distribution of Ha and 
the /-band center of the cluster (Whitmore et al. 1999). 
Such Paa shells arc diffuse and will not be identified as 
H II regions, causing the number of coincidences to be un- 
derestimated in the 5 — 7 Myr age range. For ages > 7 Myr 
the Hll region emission will not be detected as its lumi- 
nosity will be below our detection threshold (Section 5.2). 



The result of missing some (;oin(;idcnccs bcicause of the 
presence of Paa shells would be an age distribution of the 
clusters slightly younger than inferred in Section 5.3. Note 
that this effect is only relevant to the most massive clus- 
ters in Arp 299 and NGC 3256 where we expect to see 
coincidences with the present detection threshold. 

6. DISCUSSION AND CONCLUSIONS 

In this paper we have presented HST/NICMOS broad- 
band and narrow-band Paa imaging of a sample of 8 
LIRGs. The sample galaxies exhibit a range of infrared lu- 
minosities (logLiR = 10.94— 11.82 L0), as well as a variety 
of dynamical stages: isolated galaxies interacting galaxies 
and mergers. The Paa images have allowed us to identify 
the location of Hii regions, whereas the if-band contin- 
uum images have revealed the star clusters. In all galaxies 
in our sample except NGC 6240 and Zw 049.057 we have 
detected a large number of star clusters and Hii regions. 
The absolute i?-band magnitudes of the identified SSCs 
range up to Mh = —17.2 mag. A large fraction of the Hii 
region population (excluding the nuclear emission) shows 
luminosities in excess of that of 30 Doradus, the prototyp- 
ical giant H 11 region. 

The main characteristic of the Paa emission in the iso- 
lated LIRGs in our sample is the lack of strong nuclear 
emission. Most of the Hii regions are distributed along 
the spiral arms of the galaxies. The interacting/merging 
LIRGs on the other hand, show bright nuclear Paa emis- 
sion together with widely spread star formation along the 
the spiral arms and at the interface of interacting galaxies. 
The fraction of nuclear Paa emission to the total emission 
varies from system to system, as it depends on factors such 
as the age of the interaction process, the initial gas content 
of the galaxies and the relative masses of the galaxies. 

We have analyzed the properties - luminosities, sizes 
and luminosity functions of H 11 regions in LIRGs at an 
unprecedented spatial resolution (between 15 and 78 pc). 
Giant H 11 regions are ubiquitous in LIRGs and are located 
not only near the nuclei of interacting galaxies, but also 
at the interface of interacting galaxies and along the spi- 
ral arms. This population of highly luminous Hii regions 
is not observed in normal galaxies. We have fitted power 
laws to the H 11 region LFs of Arp 299 and NGC 3256 and 
found that the indices are within the values previously 
measured in the disks of normal galaxies. 

We have compared the properties of the Hii regions in 
LIRGs with a small sample of normal galaxies observed 
with the same spatial resolution and found that giant H 11 
regions are more common in LIRGs than in normal galax- 
ies. The measured sizes of giant Hii rc;gions in LIRGs 
compared to those in normal galaxies rule out the possi- 
bility that the giant H 11 regions are just aggregations of 
"normal" Hii regions. Another possibility would be that 
all Hii regions in LIRGs were extremely young. However, 
previously published spectroscopy of a few giant Hll re- 
gions in Arp 299 and NGC 3256 has shown a range of 
equivalent widths of Ha, in other words, a range of ages. 
A more plausible explanation for the presence of extremely 
luminous H II regions in LIRGs is that the regions of high 
gas pressure and density present in LIRGs, ULIRGs and 
interacting galaxies provide the necessary conditions for 
the formation of a large number of massive star (ionizing) 
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clusters. Such extreme conditions are not likely to occur 
in normal isolated galaxies. 

Despite the large numbers of near-infrared SSCs and 
H II regions identified in LIRGs, there is only a small frac- 
tion of coincidences - between 4% and 30% of the total 
number of detected H ii regions and star clusters. This is 
suggestive of an evolutionary sequence similar to that seen 
in obscured Galactic giant Hii regions. For the first few 
million years of the cluster evolution, there will be signifi- 
cant amounts of dust which may preclude the detection of 
the ionizing cluster. As the cluster ages, stellar winds and 
supernovae will dissipate part of the gas and dust, and the 
star cluster will become visible, whereas at the same time 
the H II region emission will be less luminous. 

We have used evolutionary synthesis models to repro- 
duce the observed relative fractions of young and inter- 
mediate Hii regions/clusters and old clusters in Arp 299 
and NGC 3256. Based on the observed fractions we have 
concluded that most likely the star formation occurs in in- 
stantaneous bursts rather than more extended periods of 
star formation. Using instantaneous star formation and a 
Salpeter IMF we have derived photometric masses of the 
detected star clusters of between 5 x 10'' Mq and 10® Mq. 

The fact that the peak of the ff-band luminosity oc- 
curs after approximately 9 Myr whereas at the same time 
the number of ionizing photons has dropped by about two 
orders of magnitudes from the maximum explains the lim- 
ited number of coincidences. Within the present detec- 
tion limits in Arp 299 and NGC 3256, we can only de- 
tect both Hii region emission and a star cluster for the 
most massive clusters (~ 10® Mq) during the first 7Myr. 
The near-infrared clusters with no detected Hii region 
emission will be older than approximately 7 Myr. The 
H II regions with no detected cluster counterpart are most 
likely younger than 5 Myr, and have intermediate mass 
~ 5x lO^-lO^Mo) ionizing clusters. If, as observed in ob- 
scured Galactic Hii regions, there are significant amounts 
of extinction during the first million years of the evolu- 
tion of clusters and H ii regions, then we may be missing 
the youngest star forming regions, and hence the observed 
fractions of H ii regions and coincidences will be lower lim- 
its. 

An estimate of the age distribution of the observed 
clusters can be inferred from the relative numbers of H ii 



regions, and near-infrared star clusters: the higher the 
fraction of near-infrared clusters, the older the ages of 
the detected star clusters will be. We find that the ages 
of the detected star clusters in Arp 299 and NGC 3256 
range up to 20 — 40 Myr. Older clusters possibly created 
in this or previous episodes of star formation are likely 
to exist in these systems but cannot be identified with 
the present detection threshold. Another possibility to ex- 
plain the apparent youth of the clusters in Arp 299 and 
NGC 3256 would be destruction of clusters. In that case, 
if the clusters have been created at a constant rate for the 
last ~ 100 Myr, then roughly 50% of the clusters are de- 
stroyed during that time to account for the observed frac- 
tion of clusters in these two systems. The data presented 
in this paper does not allow us to distinguish between these 
two possibilities. 

From the present observations and modelling it is clear 
that a large population of the youngest clusters (that 
is the ionizing clusters of the Hii regions with ages of 
< 5 — 6 Myr) will not be detected from near-infrared con- 
tinuum imaging alone, as only some 8 — 16% of these Hii 
regions in our sample of LIRGs (both isolated and interact- 
ing/merging systems) appear to have near-infrared cluster 
counterparts. This suggests that studies of the young star 
clusters in galaxies performed using only near-infrared con- 
tinuum imaging may be biased against the youngest star 
forming regions, which may account for up to 40% of the 
detected population of star clusters and H ii regions in 
LIRGs. 
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